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We analyze a mean-field model of electrons with pure forward scattering interac- 
tions on a square lattice which exhibits spontaneous Fermi surface symmetry break- 
ing with a d-wave order parameter: the surface expands along the k x -axis and shrinks 
along the fey-axis (or vice versa). The symmetry-broken phase is stabilized below a 
dome-shaped transition line T c ((i) , with a maximal T c near van Hove filling. The 
phase transition is usually first order at the edges of the transition line, and always 
second order around its center. The d-wave compressibility of the Fermi surface is 



however strongly enhanced even near the first order transition down to zero temper- 
ature. In the weak coupling limit the phase diagram is fully determined by a single 
non-universal energy scale, and hence dimensionless ratios of different characteristic 

in ; 

£^ , quantities are universal. Adding a uniform repulsion to the forward scattering inter- 

action, the two tricritical points at the ends of the second order transition line are 
shifted to lower temperatures. For a particularly favorable choice of hopping and 
interaction parameters one of the first order edges is replaced completely by a second 

O 

order transition line, leading to a quantum critical point. 
X ' PACS numbers: 71.18.+y, 71.10.Fd 



I. INTRODUCTION 

The low energy properties of an interacting electron system are strongly influenced by 
the shape of its Fermi surface. Residual interactions between quasi-particles near the Fermi 
surface can give rise to charge- or spin-density waves, superconductivity, or other low energy 
instabilities. Usually the Fermi surface respects the point-group symmetry of the underlying 
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lattice structure. In principle, however, electron-electron interactions can drive a Fermi 
surface deformation that breaks the orientational symmetry of the system. 

Recently, the possibility of symmetry-breaking Fermi surface deformations with a ci-wave 
order parameter, where the surface expands along the fc^-axis and shrinks along the k y -&xis 
(or vice versa), was discussed for various two-dimensional electron models on a square lattice: 
t- J,— Hubbard,^*^ and extended Hubbard^ model. The instability is driven by interactions 
in the forward scattering channel, mainly between electrons close to the van Hove points 
in the two-dimensional Brillouin zone. Referring to Pomeranchuk's& stability condition for 
isotropic Fermi liquids it has been termed "Pomeranchuk instability" by some authors. Only 
a discrete lattice symmetry is broken by the (i-wave deformation on the square lattice, in 
constrast to spontaneous Fermi surface symmetry breaking in isotropic Fermi liquids, where 
Goldstone modes play an important role. 7 The <i-wave Fermi surface deformation leads to a 
state with the same reduced symmetry as the "nematic" electron liquid defined by Kivelson 
et a/.— in a discussion of possible analogies between doped Mott insulators with charge stripe 
correlations and liquid crystal phases. 9 

Spontaneous Fermi surface symmetry breaking competes with other instabilities. In the 
slave-boson mean-field theory of the t-J model the ci-wave Fermi surface deformation is over- 
whelmed by ci-wave singlet pairing, but strongly enhanced correlations in the d-wave forward 
scattering channel and a corresponding large response to external anisotropic perturbations 
remain? 1 ^ The latter can be related to distinctive properties of magnetic excitations in 
different cuprate superconductors.—*^ In the Hubbard model near van Hove filling coexis- 
tence of superconductivity and c/-wave Fermi surface symmetry breaking has been found in 
a renormalized weak coupling perturbation expansion. 4 

Quantum critical fluctuations of the " soft" Fermi surface near a continuous zero temper- 
ature phase transition with Fermi surface symmetry breaking provide a route to non-Fermi 
liquid behavior.— Anomalously large and anisotropic quasi-particle decay rates have been 
derived for a phenomenological model, where electrons moving on a square lattice interact 
only via almost forward scattering interactions, that is only very small momentum transfers 
are allowed.— We refer to this model as the "/-model" in the following. Recently it was 
shown that the putative quantum critical point in the f-model is actually preempted by a 
first order transition, at least within mean-field theory and for various concrete choices of 
the model parameters.-^ The transition remains first order at low finite temperatures but 
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turns to second order at temperatures above a tricriticial point.— 

This paper is dedicated to a comprehensive mean-field analysis of spontaneous Fermi 
surface symmetry breaking in the f-model on a square lattice. We present results for the 
phase diagram, order parameter and Fermi surface as obtained from a numerical solution 
of the mean-field equations for various typical choices of parameters. The weak coupling 
limit is analyzed analytically. We also compute the coefficients of the Landau expansion 
of the grand canonical potential up to quartic order in the order parameter and show that 
the first order transition at low temperature is a rather robust consequence of the van Hove 
singularity in the density of states. Besides confirming the conclusions by Kee et alJ& and 
Khavkine et a/.— and providing additional numerical data especially at finite temperatures, 
we present several new results and insights. In particular, we show that in the weak coupling 
limit Fermi surface symmetry breaking is characterized by a single energy scale, which leads 
to universal behavior in terms of suitably rescaled parameters. Furthermore, we show that 
the tricritical points can be suppressed to lower temperatures by a uniform repulsion added 
to the original f-model, which, for a particularly favorable but not unphysical choice of 
hopping and interaction parameters, can even lead to a quantum critical point. Finally, 
we find that the ci-wave compressibility of the Fermi surface is usually strongly enhanced 
along the transition line down to zero temperature even if the transition is first order, which 
implies that the Fermi surface is already very soft at the transition and fluctuations should 
be important. 

The article is structured as follows. In Sec. II we define the f-model and outline the mean- 
field theory of spontaneous Fermi surface symmetry breaking. Results from a numerical 
solution of the mean-field equations are presented in Sec. III. The numerical results are 
complemented by an analysis of the Landau free energy expansion in Sec. IV and an analytic 
derivation of universal properties of the phase transition at weak coupling in Sec. V. We 
finally conclude in Sec. VI. 
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II. MODEL AND MEAN-FIELD THEORY 

We analyze the f-model on a square lattice with pure forward scattering interactions. The 
Hamiltonian reads 

H = e k ™k + ^7 ^ kk ' nk ™ k ' ^ 

k k,k' 

in standard second quantized notation, where = Y2 a n ka counts the spin-summed number 
of electrons with momentum k, and L is the number of lattice sites. For hopping amplitudes 
t, t', and t" between nearest, next-nearest, and third-nearest neighbors on the square lattice, 
respectively, the bare dispersion relation is given by 

e k = —2 [t(cos k x + cos k y ) + 2t' cos k x cos k y + i"(cos 2k x + cos 2k y )] . (2) 

The forward scattering interaction has the form 

/kk' =u — g cWk' , (3) 

with coupling constants u > and g > 0, and a function with (i^.^-wave symmetry 
such as <ik = cos k x — cos k y . This ansatz mimics the structure of the effective interaction in 
the forward scattering channel as obtained for the t- Jji Hubbard,— and extended Hubbard^ 
model. The uniform term originates directly from the repulsion between electrons and sup- 
presses the (uniform) electronic compressibility of the system. The <i-wave term enhances 
the (i-wave compressibility of the Fermi surface and drives spontaneous Fermi surface sym- 
metry breaking. In the Hubbard model it is generated by (1-loop) fluctuations, while in the 
t-J and extended Hubbard model the nearest neighbor interaction contributes directly to 
a rf-wave attraction in the forward scattering channel. For u — 0, the above model is the 
pure forward scattering limit of the f-model with small momentum transfers introduced in 

n nn 

12L Within mean-field theory, it is also equivalent to the model analyzed in Refs. 1131)141 



Ref. 



since the off-diagonal components of the quadrupole density introduced there do not affect 
the results. 

Inserting = (n^) + Sn^ into the interacting part of the model and neglecting terms of 
order (5nk) 2 , we obtain the mean-field Hamiltonian 

Hmf = y~] £kn k - ^ ^ fek (w k ) , (4) 

k k 
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where e k = e k + Se^ is a renormalized dispersion relation, which is shifted with respect to 
the bare dispersion by 

Se k = -jr^/kk' (n v ) = un + r]d k . (5) 

k' 

Here n = L^ 1 X^k( n O * s the avera g e particle density, and 

V = -j^2d k (n k ) (6) 

k 

is our order parameter, which parametrizes the amount of symmetry breaking. Note that r\ 
is real and has the dimension of energy. It vanishes as long as the momentum distribution 
function (rik) respects the symmetry of the square lattice. The grand canonical potential 
per lattice site uj = L~ l Vt is obtained from the mean-field Hamiltonian as 

-=Yg-r 2 -T EMl + e-(«^). (7) 

The stationarity conditions |^ = and |^ = (at fixed fj) yield the self-consistency 
equation for the order parameter 

k 

and the equation determining the density 

n = lE/( e k-^)> ( 9 ) 
k 

respectively, where /(£) = (e^ T + is the Fermi function. These equations follow also 
directly from the relation (n k ) = 2f(e k — fi). Note that the Eqs. (JBJ) and © are coupled for 
since n enters the self-consistency equation for 77 via e k . 
In the thermodynamic limit the mean-field theory solves the reduced version (no mo- 
mentum transfers) of the f-model exactly. One simple way to see this is by considering the 
Feynman diagrams representing the perturbation expansion of the system. All self-energy 
diagrams except the Hartree term involve integrals over momentum transfers and are thus 
suppressed at least by a factor L^ 1 . 

III. NUMERICAL RESULTS 

We first take band parameters t'/t = —1/6 and t" = 0, and solve the self-consistency 
equations eqs. (JBJ) and Q numerically. For this choice of hopping amplitudes the bare 
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dispersion relation has saddle points at (vr,0) and (0, 7r), leading to a van Hove singularity 
in the bare density of states at e° H = 4t' = —2t/3. Typical features of spontaneous Fermi 
surface symmetry breaking are captured with these parameters and are presented in the first 
three subsections. In the last subsection we investigate another set of hopping parameters, 
for which the saddle points are slightly shifted from (tt, 0) and (0, 7r), and a quantum critical 
point can be realized for suitable choices of the interaction parameters. 



A. Typical phase diagram 

We first focus on the d-wave term in Eq. © and set u = 0. A fi-T phase diagram is shown 
in Fig. ^a). The solid line denotes a second order phase transition, which turns to a first 
order transition at low T (open circles). The end points of the second order transition are 
tricritical points (solid circles), where the quadratic and quartic coefficients of the Landau 
energy expansion (see Sec. IV) vanish simultaneously. The dashed line denotes the fictitious 
second order transition that is preempted by the first order transition. The complete tran- 
sition line, T c (/i), has a dome shape. Since the Fermi surface symmetry breaking is driven 
by forward scattering of electrons mainly on the (original) Fermi surface close to the van 
Hove points,"^ the maximal T c (/i) appears for [i around the van Hove energy e° H = —2t/3; 
a slight deviation from e° H is due to finite T effects. The transition line is almost symmetric 
with respect to the fi = e° H axis. This symmetry becomes exact when t' is set to zero be- 
cause of particle-hole symmetry. Below T c the Fermi surface expands along the fc^-direction 
and shrinks along the fc y -direction or vice versa. We show results for the Fermi surface at 
low T in Figs. W[b) and (c) together with the corresponding bare Fermi surface (g = 0) for 
comparison. The Fermi surface has typically open topology in the symmetry-broken phase, 
except close to the second order transition. 

In Fig. ^d), we plot the order parameter \t]\ as a function of fi at T = 0.01£ and 0.151 
While we see a continuous transition at T = 0.15t, \rj\ exhibits a jump at T = O.Olt, 
characteristic of a first order transition. In Fig. ^e), the density n is plotted as a function of 
fi at T = 0.0 It and 0.1 5t together with that for g = 0. The density increases monotonuously 
with fi. This behavior is due to the stability condition of the system that the grand canonical 
potential must be a concave function as a function of [i, which yields an inequality, — prf > 0, 
or I 2 > 0. The density changes discontinuously at the first order transitions. The directions 
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FIG. 1: The mean-field solution for t'/t = —1/6, t"/t = 0, g/t = 1, and u = 0. (a) /i-T phase 
diagram; the total transition line, T c (fi), contains a second order transition line, T^ nd , at high 
temperatures and two first order lines, T c lst , at low temperatures; the solid circles are tricritical 
points; the dashed line, "T^ nd ", denotes a fictitious second order transition that is preempted by 
the first order transition; the dotted line indicates the van Hove energy, /i = e^ H = —2t/3. (b) and 
(c) Fermi surface in the symmetry-broken phase near the first order transition; the Fermi surface 
for g = is also shown by a gray line, (d) \x dependence of \rj\ at T = O.Olt and 0.15i. (e) \x 
dependence of n at T = O.Oli and 0.15f; the results for g = are plotted also, (f) n-T phase 
diagram; T^ nd is a second order transition temperature and solid circles are tricritical points. In 
the shaded regions, which are surrounded by T C PS , the system undergoes phase separation. 

of the density jumps are generic features required by the concavity of the grand canonical 
potential. 

In Fig. [IJf ) t ne phase diagram is plotted in the n-T plane. The second order transition 
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FIG. 2: n-T phase diagram for t'/t = -1/6, t"/t = 0, and #/t = 0.5 with (a) u = 0, (b) u/i = 1.0, 
and (c) u/t = 10. (d) u-dependence of T tr ; at higher n (solid line) and at lower n (dotted line). 

line at high T terminates at two tricritical points, below which phase separated regions, 
shown by shades, appear; dashed lines correspond to the ficticious second order transition 
shown already in Fig. ^a). Each shaded region has two phases, with different densities n\ 
and n2 (> ni). The difference between n\ and n 2 corresponds to the magnitude of the density 
jump at the first order transition point [see Fig. ^e)]. For n\ < n < n 2 the volume fraction 
of the low density phase and the high density phase is " 2 ~ n and n ~ Ul , respectively. The 

J r ° J r n2 —n\ ri2—ni ' 1 J 

former has a symmetric (symmetry-broken) Fermi surface and the latter a symmetry-broken 
(symmetric) Fermi surface in the left (right) phase separated region. 



B. Effects of uniform repulsion 

Now we switch on the uniform term in the interaction, Eq. (jSJ). The n-T phase diagrams 
are shown in Figs. Hfa)-(c) for several choices of u at g = 0.5t. Note that typical features of 
the phase diagrams are the same as in Fig. ^ although we take smaller g and finite u here. 
The second order transition line is not affected by u in the n-T phase diagram, while it 
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would be affected if plotted in the \i-T plane, because the chemical potential corresponding 
to a given density varies with u. The tricritical points, that is the end points of the second 
order transition, extend to lower T with u in favor of a second order transition, which is 
accompanied by a pronounced suppression of the width of the phase separated regions, since 
n 2 — ni is strongly reduced by the uniform term in the interaction. However, we see that 
the suppression of T tri saturates at large u, as shown in Fig. Efd), and is not strong enough 
to establish a quantum critical point. 



C. d-wave compressibility near the first order transition 

When the symmetry-broken phase is realized through a first order transition at low T, 
order parameter fluctuations are not critical at the transition. However, we now show that 
the anisotropic compressibility with a ci-wave form factor is strongly enhanced by interactions 
at the transition line, such that fluctuations can be expected to be important in spite of the 
first order character of the transition. 

The ci-wave compressibility 

dn d 

K d = -i— (10) 
dfi d 

describes the linear response of the expectation value = L _1 ^ k ci^ (n^) to the symmetry- 
breaking perturbation Hd = —fidJ2k^ n ^- The perturbation Ha induces a ci-wave shaped 
deformation of the Fermi surface. Note that the order parameter rj is directly proportional 
to rid, namely rj = —gn^. Without the form factors cik the above expressions would yield the 
conventional compressibility k = j*. For our mean field model, the ci-wave compressibility 
is given exactly by the RPA expression 

1 - gN 2 

where N 2 is a weighted density of states with d\ as a weight factor [see Eq. (fTBj) in Sec. IV] . 
The denominator 

S- l = l-gN 2 (12) 

is the inverse " Stoner factor" , which is a dimensionless measure for the enhancement of the 
ci-wave compressibility by interactions, and hence for the enhancement of order parameter 
fluctuations. 
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FIG. 3: The inverse of the Stoner enhancement factor for several choices of g and u along a first 
order transition line as sketched in the inset with the arrow; temperature is scaled by the tricritical 
temperature. 

We calculate S' 1 along one of the two first order transition lines as sketched in the inset 
of Fig. El similar results are obtained along the other side of the first order transition. The 
main panel of Fig. El shows that S^ 1 tends to zero at the tricritical temperature, that is 
the compressibility Kd diverges as expected, indicating truely critical fluctuations. At lower 
temperatures, S 1-1 is finite on the transition line. However, its value is still much smaller 
than one, especially for a smaller g; the introduction of u reinforces this tendency. The 
(i-wave compressibility is thus strongly enhanced by interactions at the first order transition 
line down to the lowest temperatures, for example by a factor of about 25 for g = 0.5t and 
u = lOt. Hence, Fermi surface and thus order parameter fluctuations can be expected to 
be important even near the first order transition at low T. In the weak coupling analysis 
presented in Sec. V we will show that S* -1 can be arbitrarily small for small g near the first 
order transition. 



= 1.0, u/t = 
= 0.5, u/t = 
g/t = 0.5, uA = 10 

f:l -I-f, ; t"/t = 




D. Quantum critical point 

The above results for our choice of hopping parameters, t'/t = —1/6 and t"/t = 0, 
are probably rather generic in the sense that other choices lead to qualitatively the same 
phase diagrams, as long as the bare kinetic energy has saddle points at (vr,0) and (0, 7r). 
Another interesting set of parameters is t'/t = —1/6 and t"/t = 1/5. For t" > \{t + 2t') 
the bare kinetic energy has saddle points at (ir ± cos -1 a, 0) and (0, it ± cos -1 a) with 
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FIG. 4: /i-T phase diagrams for several choices of u for t'/t = —1/6, t"/t = 1/5, and g/t = 0.5: 
u/t = (a), 1.0 (b), and 2.0 (c). The dotted lines indicate the van Hove energy. The lines of 
T c lst and "T^ nd " at lower u in (b) are too close to be distinguished. The \i-T phase diagrams are 
replotted function of n in (d)-(f). 

a = (t + 2t')/4t", and local minima at (tt, 0) and (0, 7r). The energy at the saddle points 
is given by e° H = -2(t - 2a 2 t"). For t'/t = -1/6 and = 1/5 one has a = § and 

fc vH — 9 t - 

A sequence of phase diagrams for t'/t = —1/6 and t"/t = 1/5 is plotted in Fig. |3J the 
coupling g is 0.5i in all cases. The transition line T c (fj,) for u = in Fig. Efa) has a dome 
shape similar to that in Fig. QJa), with a second order transition at high T and a first 
order transition at low T. Although the Fermi surface symmetry breaking is mainly due to 
electrons close to the Fermi surface near the van Hove points, the maximum of T c (/i) largely 
deviates from the van Hove energy and the dome shows a pronounced asymmetry. This is 
related to a large asymmetry of the density of states (see Sec. IV). A striking feature of the 
present parameters is a drastic suppression of one of the tricritical points by a moderate 
uniform repulsion in the interaction, which leads to a quantum critical point for sufficiently 
large u. We show fi-T phase diagrams for u/t = 1 and 2 in Figs. Efb) and (c), respectively. 
The tricritical point at higher // is only slightly suppressed by u and the suppression saturates 
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at large u as in the case of Fig- El However, the first order line at lower /i is suppressed rapidly 
with increasing u and disappears completely in favor of a continuous phase transition down to 
the lowest temperatures beyond u/t ~ 2, that is a quantum critical point is realized. Why the 
choice of hopping parameters with saddle points deviating from (n, 0) and (0, 7r) is favorable 
for a continuous phase transition at low temperatures will be discussed in connection with 
the Landau expansion in Sec. IV. 

In Figs. 0] (d)-(f), we replot the [i-T phase diagrams as a function of n. The first order 
transitions at low T lead to phase separated regions as already seen in Fig. ^ However, on 
the lower density side of the n-T phase diagram in Fig. IHe) the phase separated region opens 
disco ntinuously at the critical temperature T tr j, where the second order line terminates. In 
that case the end point of the second order line is not a tricritical point (and T tr i thus 
actually a misnomer). The quartic term of the Landau energy is positive there, but a first 
order transition nevertheless sets in due to a local minimum at finite r? in the Landau energy 
becoming a global one below T tri . 

In Fig. we show the /^-dependence of the order parameter rj for hopping parameters 
and coupling g as in Fig. HJ at a very low temperature (T/t = 0.001) and two different 
choices for u. This plot reveals that for u/t = 2 there is a first order transition within the 
symmetry-broken phase in addition and very close to the quantum critical point shown in 
Fig. 0(c). The corresponding Landau energy [Fig. E] (b)] has two minima at finite |r/|; the 
minimum at lower \t]\ has the lowest energy only for \i between the second and first order 
transition point. For larger u the first order transition in the symmetry-broken phase moves 
further away from the second order transition until it disappears completely. We show in 
Fig. E] (c) that the order parameter for u/t = 10 is continuous everywhere except at the 
first order transition at the large [i boundary of the symmetry-broken phase. The Landau 
energy for /z close to the second order transition has only one minimum as a function of 1 77 1 
in this case [Fig. El (d)]. 

No quantum critical point appears for t'/t = —1/6 and t" jt = 1/5 if the coupling constant 
g is too small or too large. In Fig. El we plot T tr j as a function of u for g/t = 1.0. Although 
T tri at lower fi is suppressed much stronger than that at higher the suppression saturates 
at a finite temperature even for large u. Remarkably the Stoner factor S is again generally 
strongly enhanced for the present hopping parameters near the transition at higher fi, while 
S is not large at low T near the transition at lower /1 for u = 0, although T tr i is lower there. 
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FIG. 5: Left panels: Order parameter \rj\ as a function of [x for hopping parameters and g as in 
Fig. temperature T/t = 0.001, and two different uniform couplings: u/t = 2 (top) and u/t = 10 
(bottom). Right panels: Landau energy uj(t]) near the phase transitions at the lower \i side of the 
corresponding order parameter plots on the left. 
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FIG. 6: u-dependence of T tr i at higher [i (solid line) and at lower fi (dotted line) for t' jt = —1/6, 
t"/t = 1/5, and g/t = 1.0. 

Obviously a low T tri does not imply a large Stoner factor. 



IV. LANDAU EXPANSION 



To gain a broader understanding of the symmetry-breaking phase transitions and their 
dependence on model parameters we consider the Landau expansion of the grand canonical 
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potential in powers of the order parameter 

u(ri) - lo{0) = l - a 2 rj 2 + i a 4 r/ 4 + . . . (13) 

up to quartic order. The function u)(rj) is given by Eq. (J7J), with the density 71(77) deter- 
mined by Eq. Q. Note that only even powers appear due to the symmetry u{— rj) = u)(rj). 
Expanding u (77) by taking 77-derivatives (see Appendix) one obtains the coefficients 

a 2 = g- 1 - N 2 (p,,T) , (14) 

a ^-<fer ) + l+M ^ fer) [ iv-ferff, (15) 

where p = // — «n(0) and 

= (16) 
k 

is a weighted density of states averaged over an energy interval of order T around p, and A^, 
Np are first and second derivatives with respect to p. Note that a 2 depends only via p on u, 
which explains the w-independence of the second order transition lines in the (n, T) phase 
diagrams in Figs. 121 and 0] The quartic coefficient 04 does not depend on g. In Fig. 0we 
plot N (p,T), N^pjT), and N%(p,,T) for T = and T = O.Olt, for the choice of hoppings 
underlying the results in Sec. IIIA-C. At zero temperature N p (p,,T) diverges logarithmically 
for p, e° H , which leads to A^(/x, T) oc —(p — e° H ) _1 and N%(p,T) oc (p — e° H ) _2 near e vH- 
These singularities are cut off for \p — e° H | < T at finite temperature. 

The critical manifold in the space spanned by g, p, and T, on which a continuous phase 
transition may occur, is given by the condition a 2 = 0, that is gN 2 (p,T) = 1. However, the 
continuous transition can be realized only if the quartic coefficient CI4 is positive. Otherwise 
it will be preempted by a first order transition. For u = 0, one has simply = —N^(p,T), 
which is obviously negative near the van Hove singularity at temperatures T \p — e° H |, 
and positive at any T > for p = e° H . Hence, the transition is first order at temperatures 
T \p — e° H |, but it can be expected to be continuous near the maximum of T c around van 
Hove filling, as is indeed the case in all numerical results. For u > there is an additional 
positive contribution to 04, which partially compensates the negative main term. However, 
the positive term is bounded by 3 [iV^ (/2, T)] 2 /Nq(P, T) even for arbitrarily large u. For p 
near e° H it is of order 1/ log \p — e° H | smaller than N%(p, T). 
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FIG. 7: The functions N (p,,T), N' 2 (ji,T), and N%(p,,T) around the van Hove energy at T = 
[(a)-(c)] and T = O.Olt [(d)-(f)] for hopping parameters t'/t = — 1/6 and t" = 0; the energy unit is 
t. 

At moderate distance from van Hove filling a numerical evaluation reveals that for hopping 
parameters t'/t = —1/6 and t" = the positive w-term can compensate a substantial amount 
of the negative term, —N^(p,T), but never turns positive. We thus confirm that a finite 
u shifts the tricritical points in the phase diagrams to lower temperatures, but does not 
produce a continuous phase transition at T = for that choice of hopping. 

The situation is very different for our second choice, t'/t = —1/6 and t"/t = 1/5. In 
that case N p (ft,T) has a step-like increase at the lower [i side of the van Hove energy at 
low T, which is generated by the local minima in the dispersion at (tt, 0) and (0, tc). This is 
illustrated for Nq in Fig- IH1 -^2 an d N± behave similarly. For p, in the step region N' 2 (fi,T) 
becomes very large at low T, while N (jJ,T) remains bounded. In the presence of a u-term 
it is thus possible to get a positive CI4 at arbitrarily low temperatures. A sizable u helps not 
only to increase 04, but also to push local minima at finite rj in u(rf) to higher energies, such 
that a continuous transition can be obtained before such a minimum becomes global. 
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FIG. 8: Nq(h,T) around the van Hove energy at T = 5 x 10~H (a) and T = O.Oli (b) for hopping 
parameters t'/t = — 1/6 and t" = 1/5; the energy unit is t. 

V. WEAK COUPLING LIMIT 



At weak coupling spontaneous Fermi surface symmetry breaking occurs only for densities 
near van Hove filling, and the transition is dominated by states with momenta near the 
saddle points of e£. In this limit the mean-field equations can be treated to a certain extent 
analytically, and the phase transition is universal in the sense that it is fully characterized 
by a single energy scale. Several ratios of physical quantities are universal dimensionless 
numbers in the weak coupling limit. 

We focus on the case u = and assume that e£ has only two degenerate saddle points 
in = (vr, 0) and ke = (0, 7r). However, the following analysis can be easily extended to 
energy bands with saddle points in other positions. Shifting energies such that e^ A = 
and choosing suitable relative momentum variables k + and one can write e£ near the 
saddle points clS db quadratic form 

1 



2m 



(17) 



where m > is a constant which can be related to the hopping amplitudes t, t', and t". The 
variables k + and /c_ are defined such that the quadratic form has no k\ and k 2 _ terms and 
that the integration measure is (27i~)~ 2 as usual. Note that the corresponding coordinate 
axes are generally not orthogonal. The momenta are restricted to a finite region around the 
saddle points by a cutoff A, that is \k±\ < A. 

If the form factor is smooth near the van Hove point it can be taken as a constant 
with alternating sign near k^ and k B , respectively, such as 



1 for k « k,i , 



•1 for k k 



(18) 



B 
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Any other constant (7^ 1) could be absorbed in the coupling g. The self-consistency condition 
Eq. (jHJ) becomes simply 

V = 9 ( n B - n A ) , (19) 

where 

n A ,B = 2[ ^/(4±»?-A*) (20) 
Ja,b W 

is the contribution of the momentum space region around and kg (limited by A) to the 
density 



A. Ground state 

At T = the integral over the Fermi function is simply the volume of occupied states. 
An elementary integration yields 

A 2 2m / A 2 \ 

n A ,B = — + — G" =F V) f 1 + log— -log|/i=F^| J • ( 21 ) 

The self-consistency equation can thus be written as 

V = 9[(p-v) ]o s\^-'n\-(P' + V) log \fi + 77 1 + 2rj (1 + log e A ) ] (22) 

with the dimensionless coupling g = and the cutoff energy e\ = A 2 /2m. Integrating 
duj/dri = <7 x ?7 — (% — Ua) over 77 one obtains the 77-dependence of the grand canonical 
potential 



2m 



n 2 



1 1 H , \ 

1 + loge A 

2g 2 V B ' 



rf 



+ \ + V) 2 log |/i + 77I + ^ (fi - r] f log I)" - rj\ J + const - > (23) 

where the constant does not depend on 77 and can be chosen such that uj(Q) = 0. 

For i-i = 0, where the bare Fermi surface touches the van Hove points, the self-consistency 
equation for 77 becomes simply 

77 = -2(7 7? log \r]\ + 2(7 77 log(ee A ) • (24) 

Besides the trivial solution 77 = this equation has the two degenerate solutions 

77o = ±e6 A e- 1 /( 2 §). (25) 
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It is easy to see that the solution 77 = corresponds to a maximum of u>(rj) and is thus 
unstable. The total density remains unaffected by the symmetry breaking, since 5ns = Sn^ 
for /n = 0. This is also true to high accuracy for the density at van Hove filling for the 
numerical result shown in Fig. ^e). Note that at van Hove filling symmetry breaking occurs 
for arbitrarily small g, which is due to the logarithmic divergence of the density of states. 

For |i / we introduce a rescaled order parameter fj = r\j\i. The self-consistency 
condition for rj can be written in terms of fj as 

? = (1 - fj) log 1 1 - fj\ - (1 + fj) log 1 1 + fj\ (26) 

with a renormalized coupling constant given by 

I = - + 21og|/i| -2(l + loge A ) . (27) 
9 9 

The grand canonical potential can be written as uj{rf) = ^ {i 2 u(fj), where 

&{v) = Q~ + ^(l + ^) 2 log|l + r/| + ^(l-r/) 2 log|l-77| . (28) 

Note that \i and the cutoff have been completely absorbed in the renormalized coupling g. 
The rescaled self-consistency equation and u(fj) are universal in the sense that they depend 
only via g on all input parameters. 

Minimizing G){fj) one finds that a first order transition occurs at the universal critical 
coupling g~i pa —0.692 with a universal jump of the dimensionless order parameter 

= -PH w 1.720 (29) 
[Mil 

The Fermi surface thus opens immediately at the transition. Inverting the relation between 
g and g one obtains the //-dependence of g~\ as 

5lM =l + 2, l( log- + l) 2l^f' (30) 

Note that g\ > for e\ although g\ < 0. For fixed g, on the other hand, the critical 

value of fi at which the first order transition occurs is 

M=e l+l/( 2i h) eAe -l/(29). (31) 

The curvature of Cj{fj) at fj — becomes negative only for g > g 2 = —0.5. At fixed \i this 
requires couplings g > g^A*) > gi(^)- For given g, the critical value of \i for a continuous 
transition is 

\^\ = e A e- 1 ^\ (32) 
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which is smaller than \fii\. Hence, the critical point for a continuous transition is not 
reached, but preempted by a first order transition. The ratio | | / 1 /x 2 1 = e 1+1 ^ 2 ^ ~ 1.320 
is a universal number. 

In the symmetric state, the interaction induced enhancement of the <i-wave compressibility 
of the Fermi surface is given by the "Stoner factor" S — [1 — gN 2 (/i)]" 1 , which is related 
to the quadratic coefficient of the Landau expansion by S = (ga 2 )~ 1 , see Eq. (JT3J). Since 
a 2 vanishes for g = g 2 , one has a 2 = g^ 1 — g 2 ■ At the first order transition in the ground 
state, a 2 is given by 

a,\ x = QT 1 ~ 9? = ^ fe" 1 - & 1 ) « 0.555 ^ , (33) 

where we have used Eq. f!27|) in the second step. The <i-wave compressibility is thus enhanced 
by a factor 

s ' = oil (34) 

at the first order transition. For a weak coupling g <C 1, this enhancement is very large. 



B. Finite temperature 

We now compute two characteric temperature scales in the weak coupling limit, namely 
the transition temperature at van Hove filling, To, and the tricritical temperature, T tr j. To 
this end we write the functions N p (ft,T) in the form 

N p (jm, T) = - J deN p (e) f'(e - fj) , (35) 

where N p (e) = N p (e,0), which, for the quadratic dispersion in Eq. (|T7|) with a cutoff A, is 
given by 

iV p (e) = ^log^ (36) 

with |e| < 6a- Note that ft = /i for u = 0. 

The transition temperature at van Hove filling is obtained by setting the quadratic co- 
efficient in the Landau expansion a 2 to zero at /i = 0, that is by solving the equation 
gN 2 (0,T ) = 1. Using 

" de log M f(e) log | - log \ + 7 , (37) 
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where 7 ~ 0.577 is the Euler constant, one obtains 

T = — e A e-^. (38) 

7T 

Note that the numerical prefactor coincides precisely with the one in the BCS formula for 
the critical temperature of a weak coupling superconductor. 

At the tricritical point a 2 and 04 both vanish. The tricritical temperature T tr ; and the cor- 
responding chemical potential // tr i are thus determined by the two equations gA 2 (/itri, T tr i) = 
1 and A 7 "^' (/i tri , T tr i) = 0. Using Eqs. and (f3T)|) one obtains 

gN 2 (fx,T) eA ^°° 23[log^ + a(/2)] , (39) 

where ft = fi/T and the dimensionless function a(jl) is defined as 

f°° d 1 

Furthermore 

N'l^T) J^L&(/i) (41) 

with 

6(0 = j jx log + (42) 

The latter function vanishes for /i = ±/i tr i with /z t ri ~ 1.911. Setting the right hand side of 
Eq. (|3TIJ) equal to one and solving for T then yields 

T tri = e- a e A e- 1 /^ (43) 

where a = a(fl tri ) « 0.4515. Hence, the tricritical temperature and the critical temperature 
at van Hove filling form the universal ratio 

= ~ 0.5614 . (44) 

C. Comparison with numerical results 

Above we have computed several physical quantities characterizing the phase transition, 
which are all proportional to the same energy scale e\e~ 1 ^ 2 ^\ with universal prefactors in 
the weak coupling limit. Hence ratios of these quantities are universal numbers. 
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We have checked universal ratios against results from the numerical solution of the mean- 
field equations for hopping parameters t'/t = —1/6, t" = and coupling u — 0. At zero 
temperature we have checked ratios involving the order parameter t]q at van Hove filling, the 
order parameter jump r]\ at the first order transition, and the distances of \L\ and fi2 from 
e° H . For a comparison of the order parameter one has to take into account a factor two due 
to the different size of d^ near the van Hove points, which is ±2 in the numerical calculation 
with dk = cosk x — cosk y , but ±1 in the weak coupling model. At finite temperature we 
have compared the transition temperature T at van Hove filling and the tricritical point. 
For g/t = 0.5 all ratios agree within one percent error with the predicted universal numbers. 
For the stronger coupling g/t = 1 the deviation increases to two or three percent at zero 
temperature and up to around five percent at finite temperature. 



VI. CONCLUSION 



In summary, we have analyzed a mean-field model for Fermi surface symmetry breaking 
with a rf-wave order parameter on a square lattice. We have confirmed the qualitative 
properties of the phase diagram reported already by Kee et a/.— and Khavkine et alM We 
have provided further numerical evidence and analytic arguments showing that the phase 
transition is typically first order at low temperatures. This implies that a stability analysis 
of microscopic models, for example by renormalization group methods, should not focus on 
diverging susceptibilities only. At weak coupling the transition is fully characterized by a 
single energy scale and can thus be described by universal dimensionless functions of suitably 
rescaled parameters, which leads to various universal ratios of different quantities. The 
tricritical points separating first and second order behavior are shifted to lower temperatures 
by adding a repulsive constant contribution to the forward scattering interaction, which, for 
a favorable choice of hopping and interaction parameters, can even lead to a quantum critical 
point. Although the phase transition is usually first order at low T, we have found that the 
c?-wave compressibility at the transition can be enhanced significantly by interactions, which 
implies that Fermi surface fluctuations induced by interactions with a small finite momentum 
transfer are expected to be important even near the first order transition. The role of 
fluctuations, in particular their influence on the phase transition, remains an interesting 
subject for future studies. 
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APPENDIX A: DERIVATION OF LANDAU EXPANSION 



Here we derive the expressions Eqs. (114)1 and (jl5j) for the coefficients of the Landau 
expansion of Lj(n), by taking derivatives with respect to rj at fixed ]i and with n(rj) determined 
by Eq. (JHJ). The first derivative is 



duj duj ^du dn n ^ 2 y-^ ^ _ 
dr/ dn dn dr] g L 



(Al) 



where we have used the stationarity condition = 0. The second derivative is 



f)n 



dn 2 



9 



L 



dl + d k u 



dn 
dn 



f(e k - M) • 



(A2) 



Exploiting the symmetry n(n) = n(—n), which implies that odd derivatives of n(n) vanish 
at rj = 0, one obtains 



d 2 u 



a 2 



dn 2 



g- 1 + - dl f(ei + u n(0) - /x) = g~ l - N 2 (fi, T) 



(A3) 



•7=0 - k 

with p, = u — un(0). 

Differentiating twice more and setting n = 0, one obtains 



d 4 u 



a 4 



dn 4 



-N':(fx,T) + 3uN^,T) 



d 2 n 



ri=0 



drf 



ri=0 



Applying two ^-derivatives to n(r]) as given by Eq. Q we get 



d 2 n 



drf 



d 2 n 



ri=0 



N' 2 {-u,T)-uN {-u,T) 



(A4) 



(A5) 



»7=0 



Solving for n"(0) and inserting into Eq. (|A4|) we obtain Eq. f)15j) . 
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